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Abstract 



fTj ' Gamma-Ray Bursts (GRBs) are expected to efficiently accelerate protons up to 

relativistic energies. High-energy photons can originate from decay of neutral pi- 
►^ , ons produced by the interaction of these protons with the medium surrounding 

Q^ ' the source. In the same hadronic chain, high energy neutrinos are expected to be 

OC ■ produced from decay of charged pions, as well as from decay of generated muons. 

h.^^ . Neutrinos can travel cosmological distances being much less absorbed than photons, 

f^ I thus providing a powerful tool for the investigation of cosmic sources. In the frame 

of a recently proposed hadronic model for the emission of high-energy 7-rays, we es- 
timate the neutrino emission from GRBs and calculate the neutrino fluxes at Earth 
as a function of the model parameters. The detectability of the expected signal by 
current and future experiments is also extensively discussed. 



o 
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1 Introduction 



Nowadays, it is widely recognized that GRBs have cosmological origin, as sug- 
gested by Beppo-SAX observations [1] of counterparts in the X-ray band of 
some GRBs. Thanks to these observations, it has been also possible to deter- 
mine with great accuracy the GRB position and consequently, through point- 
ing of optical and radio telescopes, to evaluate the distance of some GRBs. 
The derived amount of energy radiated by assuming an isotropic explosion 
can result in some cases of the order of 6^'''° ~ 10^^ — 10^^ erg in the soft 7-ray 
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(50 - 300 keV) energy band, showing that in a GRB event up to one solar 
mass can be radiated in form of soft 7-rays. Note that if GRBs are beamed 
sources as indicated by afterglow observations [2], the emitted energy would 
be decreased by the factor Af2/47r, where AQ is the beaming angle. 

According to the fireball model for GRBs [3] , relativistic electrons are acceler- 
ated at shock fronts occurring near the interface of the expanding relativistic 
shells (external shock models), or at shocks forming within the unsteady out- 
flow itself (internal shock models). Soft 7-rays are then produced through 
leptonic processes, such as synchrotron or Inverse Compton emission. 

Many GRBs have been recognized as high-energy 7-rays emitters by the 
EGRET experiment which detected 7-rays in the energy band between 30 
MeV and 20 GeV from a high percentage of the bright bursts found by BATSE 
[4-6]. These observations clearly show that the GRB phenomenon is not ex- 
clusively the domain of soft-energy leptonic processes but may include high- 
energy hadronic mechanisms. 

It has been suggested (see e.g. [7]) that highly relativistic protons accelerated 
at shock front can produce quite efficiently very high-energy photons and neu- 
trinos through p7 interactions within the dense photon field arising from the 
GRB [8], or by mutual pp interactions in the ejecta [9], or by pn inelastic 
collisions with neutrons, if present in the original explosion [10] . 
The emission of high-energy neutrinos and 7-rays has been also investigated 
within the so-called cannonball scenario, by which jets of highly relativistic 
balls emitted in supernovae explosions collide with the surrounding shell pro- 
ducing quasi thermal radiation (due to the collision heating), and high energy 
7-rays and neutrinos (through pp interaction in the shell) [11] 

Very recently, a model for the production of high-energy 7-rays in GRBs with 
Ej > 1 GeV has been proposed [12]: high-energy protons emerging from the 
GRB source may subsequently interact with external nucleons, giving rise 
to neutral pions, which in turn promptly decay into high-energy photons. A 
key requirement of the model is the existence of a dense enough cloud (as a 
target ior pN interactions), which surrounds the GRB source, or, alternatively, 
situated along the line of sigh to the observer and near enough the source. 
A delay in the arrival times of the high-energy 7-rays with respect to the 
soft 7-rays and a temporal spread of the signal at Earth is predicted as a 
consequence of the proton defiection due to the presence of magnetic fields 
in the region surrounding the GRB source. Clearly, in addition to neutral 
pion production, inelastic pN collisions also produce charged pions, kaons 
and small quantities of other mesons. If their main free path for interaction 
Aj = (^atct^at)"^, is much larger than their decay path A^ = Ttc (with T = 
E'^/m^c^), they decay into charged leptons and neutrinos [14] \^ Here n^ is 



Note that even in the case of a cloud with density as high as n^ — 10 cm "^^ the 



the nucleoli number density, cr,rAr is the pion-nucleon cross section and r is 
the pion proper mean hfe. The description of this model is given with some 
detail in section 2. In section 3, we estimate the neutrino emission from GRBs 
as expected within this frame and compare our results with the theoretical 
predictions by other models for the production of neutrinos from astrophysical 
sources. In section 4 we finally discuss the detectability of the expected signal 
in terms of upward-going neutrino-induced muons by taking into account the 
experimental constrains given by underground and underwater/ice detectors. 



2 The Model 

According to [12], we assume that the GRB source emits an amount of energy 
£p ~ lO^'^Ail erg in form of relativistic protons [13], released at a distance 
To ~ 10^^ cm from the central source, during a time At ~ 1 s. 

Moreover, we require that there exists a dense enough cloud near/around 
the GRB source (with un ~ 10^ — 10^^ cm~^) so that accelerated protons 
subsequently interact with the nucleons in the cloud {apN ~ 3 x 10~^^ cm^), 
giving rise to neutral and charged pions and ultimately to high-energy photons 
and neutrinos 0. 

The initial proton flux Jp{Ep,r) is taken to be a function of both the proton 
energy Ep and the radial coordinate r as in the following: 

Jp{Ep,r) = C (Ep/GeV)-"^/(r) [protons cm'^ s~i sr'^ GeV~^] . (1) 

Here, the radial profile, for r > tq, has the form 



— exp[-crpAr TiAf /(Ep,r)] if I < Rl{Ep) 
f{r) = { \r J 

if / > RiiE, 
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where RiiEp) = Ep/{eB) is the proton Larmor radius, l{Ep, r) = a{Ep, r)RL{Ep 
is the path-length travelled by protons starting at tq up to their first pN in- 
teraction occurring at distance r from the source and a{Ep,r) is the proton 



condition A^/Aj << 1 is well verified, so that all the produced secondary unstable 
mesons decay before interacting with the cloud itself. 

^ The average length travelled by protons in the surrounding medium results to be 
AR = {apNiiN)^^ ~ 10^'^ pc (10"'^'' cm^^/njv). Correspondingly, the mass enclosed 
will be Mci ~ 0.4 M0(ro/lO^^ cm)^ in a shell-like cloud geometry around the GRB 
and Mci ~ 0.2 Mq{1Q^^ cva~^ /hn)'^ in the case of an intervening cloud near the 
GRB source. 



deflection angle due to the intervening magnetic fleld B (tipically B ~ l/xG) 
given by: 

(r — tq \ 



As far as the value of the energy spectral index a^ is concerned, since protons 
are accelerated by shock waves moving with ultra-relativistic velocities (corre- 
sponding to GRB Lorentz factor T ~ 300, see e. g. [15] and diffusive processes 
are negligible during the time before pN interactions, we adopt the relevant 
value ap~2.2 [16,17]. 



Finally, the constant C in eq. (1), results from the normalization condition 
where Eq = TiripC^ is the minimum proton energy. 



m rl j^ J,{E.p, ro) Ep dE, = ^^ , (3) 



Following Dermer [18], the production spectrum of high-energy 7-rays result- 
ing from pN collisions, via 7r° -^ 77 decay, is given by the source function 



f-t-OO f + OC 

q^{E^,r) = IQix'^riN / dE^, / dEp Jp{Ep,r)x 

(4) 
f' dcose E*^ [7 cm-3 s-i GeV-i] , 

JcOSdmn^ dV 



where 6 is the angle between the line of sight and the proton direction, asterisks 
refer to CM frame and E*d^a* /dp*^ is the Lorentz invariant cross-section for 
7r° production in high-energy {Ep < 10'^ GeV) pN collisions [^j. 

Numerical values of q^{Ey,r) depend on the main model parameters, namely 
the GRB energetics £p, the nucleon number density rij^ and the magnetic fleld 
strength B (for more details see [12]). 

An important effect of the presence of magnetic flelds is a temporal spread of 
the emitted high-energy 7-ray signal on Earth. To better analyze this point, 
we can consider the proton injection as instantaneous, since the acceleration 
time At ~ 1 s is much smaller than the average proton interaction time 
(cTpAf^Tvc)"^ ~ 10^ s (10"*^° cm~^/niv). Accordingly, denoting by to the instant 
at which photons arrive on Earth in the case B = (no proton deflection). 



^ Mori [19] had verified that extrapolation up to Ep ~ 10^ GeV of the invariant 
cross-section used by Dernier [18] is in reasonable agreement with Monte Carlo 
simulations by several pN event generators. 



the delay ti — to can be expressed as a function of the radial coordinate r^ at 
which pN interactions occur and the proton Larmor radius Rl{Ep): 

RiiEp) ( . Ti-ro ri-ro\ 

H —tn = ^— arcsm — — - . (5) 

c V MEp) Rl{Ep)) ^ ' 

Vice versa, one can say that photons which arrive on Earth within the time 
window {tQ.ti) are produced by protons of energy Ep only if these protons have 
interacted with nucleons in the shell between r^ and ri{Ep). 

Correspondingly, the source spectrum Q^{E^;tQ,ti) of 7-rays produced in a 
GRB event within the temporal window ti — to is given by 

/•+00 /■+00 

Q^{EY,to,ti) = IGn'^UN At/ dE,, dEp x 

JE-rp-i" J En 



/ drJp{Ep,r) exp[-apNnNl{Ep,r)] x (6) 

Jra 



[' dcose E*% [7 GeV-1] , 

with the proton injection time At ~ 1 s. 

In Fig. 1, assuming different values for the model parameters un- and B 
in the range 10^ — 10^^ cm~'^ and 10"^ — 10'^ G, respectively, we show 
Q-y{E^ ; 0, 200 s) in the first 200 s, starting from the inset of the signal 
to = 0. From Fig. 1, one can see that for a given photon energy E^ the 
differential 7-ray production increases with increasing nj^f and decreasing B. 
The effect is particularly important at low E^ values. In fact, low-energy pho- 
tons are mainly generated by low-energy protons with Ep roughly one order 
of magnitude above E^: they are more easily deflected by the intervening 
magnetic fields giving rise to photons that cannot reach the observer since 
a{Ep,r) > 7r/2. However, our numerical results show that in a large range of 
values for n^r and B, the decrease of Q^{E^;tQ,ti) for increasing values of B 
can be counterbalanced by increasing n^v- Moreover, for each assumed value 
of B, there exists a limiting density (n^r — 10^^ cm~'^ in the case 5 ~ 1 /iG) 
above which most of the protons interact before suffer substantial deflection. 

Note that in deriving eq. (6) we have neglected re-absorption of the produced 
7-rays in the cloud itself 0. Clearly, re-absorption of neutrinos within the 



'^ High-energy 7-rays may in principle be re-absorbed through the interactions with 
the atomic hydrogen, provided that the path they have to travel inside the cloud 
(after being produced) is long enough. However, clouds nearby the GRB source 
are expected to present a high ionization degree, so that the relevant absorption 
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Fig. 1. The source 7-ray spectrum Q-y{E-y;0, 200s) in the first 200 s of the signal 
life for different values of the cloud number density n^ = 10^ — 10^^ cm~^ and of the 
magnetic field strength B = 10^^ — 10"^ G. The GRB total energetics is assumed 
to be 10^^ erg/sr. 

cloud itself can be completely neglected. 



3 High-energy neutrino emission 



It is well known that, for a proton power law spectrum Jp{Ep) oc E^^-p and 
inclusive cross sections for meson production obeying Feynman scaling, the 
produced spectrum of 7-rays, u^s, z/^'s and z/t-'s are all proportional at high 
energies (see e. g. [20]) Following the calculation developed by Dar [21] and 
by Lipari [22], we assume in particular that muon neutrino and gamma ray 



spectra by pN interactions are related through Qi, 



KQ^ with K ~ 0.7. 



From eq. (6), we can then estimate the differential neutrino fiuences on Earth 



processes are je and jp interactions, which imply a large 7 free-path before re- 
absorption. 



within the time interval ti — to 

Q,[E,{l + zy,to,U] (l+z) 



Fy[Ey; tQ,ti 



AnDliz) 



[u cm-2 GeV-^] , (7) 



where z is the GRB redshift and 



Dl{z) 



Hoq'o 



zqq + (go 



l)(-l + V'2go^+l)) 



is the GRB luminosity distance, with Hq ~ 70 km s~^ Mpc^^ and go — 0.15. 
In Fig. 2, assuming z = 1 and Sp = 10^^ erg/sr, the neutrino fluence F^^E^,; 
0, 200 s) in the initial 200 s is given as a function of neutrino energy for 
selected values of njv and B parameters. A comparison with Fig. 1 shows 
that, apart the energy z shift, the neutrino fluence follows, as a function of 
model parameters, the same behaviour of the source 7-ray spectrum. 
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Fig. 2. Expected f^ + i?f^ neutrino fluences Fy{Ej^; 0, 200 s) in the first 200 s of the 
signal life for different values of the cloud number density tin = 10^ — 10^^ cm~^ 
and of the magnetic field strength B = 10"^ — 10"'^ G. The GRB total energetics 
and redshift are assumed to be 10^^ erg/sr and z = 1, respectively. 



Starting from the neutrino fluence F^{E^; to, ti) during the time interval ti — to, 
one can estimate an average neutrino flux on Earth from a single source as 



follows: 



ti — to 



Otherwise, it is also possible that neutrinos from unresolved GRBs sources 
give a contribution to the diffuse high energy neutrino background. In this re- 
spect, due to the widely accepted cosmological origin of GRBs, we adopted, as 
recently done in [23] , the cosmological evolution scheme developed by Schmidt 
[24] . In particular, we assume that the rate of ORB events, per unit comoving 
volume and per unit comoving time, is given by 

ho{l + zy ioT < z < Zp 

n[z] = < 

no(l + ZpY for z > Zp 

with Zp = 1 and p = 3.32. The parameter tiq is then determined so that the 
observed (in the soft 7-ray band) total GRB event rate 

Ngrb = / T^r^ir ^-2 10 

JO 1 + z dz 

reads Ngrb — 10'^ yr^^. Here, dV = AnDl^dDM is the comoving volume 
element. Dm = Dl/(1 + z) is the proper motion distance, related to the 
luminosity distance D^ given in eq. (8), z^ax is taken equal 5. Moreover, in 
eq. (10) we also take into account that the time interval between bursts is 
stretched by the factor {1 + z). 

For simplicity, we assume that GRBs are standard candles, each source with 
£p = 10^^ erg/sr, B = l/iG and n^ = 10^ cm~^. Correspondingly, the diffuse 
neutrino flux expected on Earth results to be 

^ p Q„K(14-.)0.+oo] n(^,V ^^ 

Jo 4:nDf^{z) 1 + z dz 

In Fig. 3, ^,y{E,,) is given as a function of Ei^. For comparison we show the 
diffuse neutrino flux from GRB events (via p~f interactions) predicted by the 
Waxman & Bahcall model [8] as well as the upper bound implied by high 
energy cosmic ray observations estimated by the same authors [25] . Moreover, 
in Fig. 3, we give the atmospheric neutrino flux averaged on the zenith angle as 
calculated by Gaisser et al. [26]. Actually, at energy above 10 TeV, the isotropic 
contribution of neutrinos produced in the atmosphere by the semileptonic 
decays of charmed hadrons (prompt neutrinos), starts to dominate over the 
conventional pion and kaon decay induced atmospheric neutrino flux. In Fig. 3 
we report the contribution of prompt neutrinos as calculated in [27] within the 
frame of the Recombination Quark Parton Model. It should be remarked that 
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Fig. 3. Diffuse fluxes of u^ + U^ from this model assuming no evolution [GRB pN 
(1)] and Schmidt's evolution [GRB pN (2)], respectively. The predictions according 
to other models of AGNs and GRBs are also shown: SDSS [31], W [25]. The dashed 
curve refers to the angle average atmospheric neutrino flux [26] ; the dash-dotted line 
refers to the flux of prompt neutrinos according to one of the highest predictions 
given in [27]; the dotted line is the theoretical upper bound to neutrino flux from 
astrophysical sources as calculated in [25]. Some current experimental upper limits 
are also shown: AMANDA [32], MACRO [33], Frejus [28] and BAIKAL [30]. 



there is some theoretical uncertainty on the calculation of prompt neutrinos 
flux, as well as on the energy at which this contribution crosses over the 
conventional atmospheric flux. In Fig. 3, we considered one of the highest 
prediction among those given in [27]. Finally, the upper limits from current 
experiments (Frejus [28], AMANDA [29] and BAIKAL [30]) are shown, derived 
by assuming a power law spectrum (oc -E~^) an initial neutrino flux. 



4 Neutrino signal detectability 



4.1 Neutrino propagation 



High-energy neutrinos from astrophysical sources can be detected by collecting 
the upward-going neutrino-induced muons produced through charge-current 
interactions with the matter surrounding the detector. Based upon this de- 
tection principle, current underground and underwater/ice experiments are 
surveying the sky in order to search for neutrino sources. In the frame of the 
model described above, we have investigated by numerical calculation the de- 
tectability of the signal expected either from a single point-like GRB, or from 
the diffuse flux due to unresolved GRBs. Neutrino propagation through the 
Earth has been developed following the approach discussed in [34] , where they 
give an analytical method for the precise calculation of the spectrum of high 
energy neutrinos propagating through a dense medium of any thickness, for 
initial spectrum and for cross sections of any form. 

As a preliminary observation, it should be noticed that the results of atmo- 
spheric neutrino experiments [35,36] suggest that muon neutrinos can oscillate 
into tau neutrinos, thus affecting the expected muon event rate on the detector. 
This effect has been extensively discussed e.g. in [37]. Here, flavour oscillation 
will be not considered. 

At energy above 1 GeV neutrino-matter interaction is dominated by the pro- 
cess of deep inelastic scattering (DIS) on nucleons (charge and neutral current 
processes), since the contributions of both elastic and quasi-elastic interactions 
become negligible. The main difficulty for the evaluation of the deep inelastic 
neutrino- nucleon cross sections arises from the fact that, in spite of the lack of 
experimental data, at very high energy the details of nucleon structure become 
important. 

Detectable muons are produced through charge current interaction as in eq. (12). 
Neutral current interaction described by eq. (13) causes instead a modulation 
in the spectrum of the interacting neutrinos. 

u^{u^) + N-.fi-{fi+)+X , (12) 

u^(u^) + N ^ u^iu^) + X' . (13) 

The total cross section for these processes can be written in terms of differential 
cross sections as follows: 

1 _!!!£. J cc 

<otiE.) = l'""^iK,y)dy, (14) 

<t{E.)=l -j^{E^,y)dy, (15) 

where E,^ is the energy of the incoming neutrino, m^ is the muon mass and 
y is defined by Ei/E^, with Ei = E^ — E^ for charge current interaction [E^ 



10 



is the energy of the outcoming muon) and Ei = E^ — E^ for neutral current 
interaction [E^ is the energy of the outcoming neutrino). 
Let Fy{Ey,x) be the differential energy spectrum of neutrinos at a column 
depth X = /q p{L')dL', where p{L) is the density of the medium at dis- 
tance L from the boundary, measured along the neutrino beam path. The 
one-dimensional transport equation can be expressed in the form [34]: 



dx \i,{Ey) Jo An dy 

with the boundary condition F^{E^,0) = F^{Ey). Here, \u{E^) is the the 
neutrino interaction length in g cm~^ defined by: 

^''^^''^ ^ {NA/AN){aUE,) + aUE,)) ' ^^^^ 

where N^ is the Avogadro number and A^ = 1 g mol~^. 

The first term in the right side of eq. (16) accounts for the absorption of 
neutrinos of energy Ey^ due to charge and neutral current interactions; the 
second term describes the additive contribution of higher energy neutrinos 
"regenerating" through neutral current interaction at energy E^ . The solution 
of eq. (16) is taken of the form 



F,(E,,a:)=F°(E,)exp 



X 



MEu,x) 



where the function Ai,{Ei^,x) is the effective absorbtion length containing the 
relevant information on the modulation of the neutrino flux Fj^{Ey^x) along 
the propagation (see Ref. [34] for further details). The differential v^j^N and 
V^N cross sections have been calculated by using the approach in Ref. [38] 
based on the renormalizat ion-group- improved parton model and on the last 
available experimental information about the quark structure of the nucleon. 
They are written in terms of the Bjorken scaling variables y and x = Q'^/2MEi, 
where —Q^ is the invariant momentum transferred between the incoming neu- 
trino and the outgoing lepton (muon or neutrino, for charge and neutral cur- 
rent interaction, respectively). The detailed formulas used for the calculation 
are given in [38]. 

Various versions of different sets of parton density functions q{x, Q"^) are col- 
lected in a large CERN software library PDFLIB [39]; they can be simply 
accessed by setting few parameters to choose the desired version. In this calcu- 
lations, we selected the third version of the CTEQ collaboration model [40] for 
the next-to-leading order (NLO) quark distributions in the deep-inelastic scat- 
tering factorization scheme. The CTEQ3 distributions are particularly suit- 
able for high energy calculations since the numerical evolution is provided for 
very low x. Unfortunately, the uncertainty due to this extrapolation may be 
rather large and is hard to estimate. The Q^ evolution is realized by the NLO 
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Gribov-Lipatov-Altarelli-Parisi equations from initial Qg — 2-56 GeV^. (The 
new version of parton density functions by CTEQ group [41] doesn't produce 
significant changes for this calculation). 

Fig. 4 shows the neutrino interaction length \ij{Ey) for deep inelastic scat- 
tering off an isoscalar nucleon A^ = {n + p)/2 (n ^ neutron, p -^ proton), 
for v^ (solid line) and v^ (dashed line). The dotted line placed at the level of 
the Earth diameter demonstrates that the Earth becomes opaque to neutrinos 
with energy exceeding the TeV energy-range. 
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Fig. 4. Neutrino interaction length for deep inelastic scattering off an isoscalar nu- 
cleon N = {n + p)/2, for u^ (solid line) and v^ (dashed line). The calculation has 
been done for the CTEQ3 set of parton density functions [40] . The dotted line gives 
the Earth diameter. 



-^.^ Muon propagation 



High energy muons propagating through matter mainly interact by quasi- 
continuous (ionization) and discrete {e^e~ pair production, bremsstrahlung 
and photonuclear interaction) electromagnetic scattering. Ionization domi- 
nates at energy lower than few hundreds of GeV , whereas above the TeV 
region energy losses are mainly due to radiative processes. 
At this stage, muon propagation has been treated analytically, according to 
the approach suggested in [42]. Following [42], the one- dimensional transport 
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equation can be written as: 

"^ -{f3{E^)F,{E^,x))+g{E„x), (19) 



dF^iE^,x) _ d 



with the boundary condition 

F,{E„0) = 0. 

Here, F^{E^^ x) is the muon flux at energy E^ and depth x, P{E^) is the aver- 
age total energy loss per unit length due to electromagnetic muon scattering 
(ionization plus radiative processes) and ^(-E^, x) is the "source" function in- 
troduced to describe the production of muons through neutrino charge current 
interaction. Q{Efj_, x) is given by: 



An Jo dy \l-y J 



\^.A ^. (20) 

The solution of eq. (19) can be written in the form: 

F,{E„x) = J^^ .^^g{E^,x-n{E^,E,)), (21) 

where the function TZ{E^, E^) defined by 



'^(^•^'■'-rji) '^^' 



represents the mean depth crossed by a muon of initial energy E and final 
energy E^; £{E^, x) labels the initial energy needed by a muon to reach the 
depth X with final energy E^^. The calculation of muon energy losses is crucial 
to correctly evaluate the expected muon fiux on the detector. The muon energy 
losses due to radiative processes (e.g. e+e~ pair production, bremsstrahlung 
and photonuclear interaction) can be expressed in terms of differential cross 
section for the specific muonic interaction process k as follows: 



dE\ Na „ r^'^^ da 



/dE\ 



k 



^M / v--{v,E^)dv, (23) 



\ dx / k A Jvjnin dV 

where Na and A are Avogadro's number and the mass number respectively; v 
is the fraction of initial energy E^ lost by muon at the occurrence of the process 
k] Vmin and Vmax are the kinematical limits for the allowed values of v; x is 
the thickness of crossed matter, expressed in g cm~^. The total average muon 
energy losses I3{E^) results from the sum of the radiative interaction terms 
plus the quasi-continuous contribution due to ionization. The explicit formulas 
used in this work are given in [43], accordingly to the theoretical framework 
discussed in the standard reference by Lohmann and Voss [44]. Fig. 5 shows 
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the total average muon energy loss (thick line) and the contributions from the 
single electromagnetic process. 




Log,o'(E (GeV))' 



Fig. 5. Total average muon energy losses in standard rock (thick line). The contri- 
butions of each electromagnetic process have been also plotted, ion -^ ionization, 
pp — > e^e~ pair production, 6s -^ bremsstrahlung, ph — > photonuclear interaction. 



4.3 Discussion and results 



In Fig. 6, assuming the same parameter values adopted in Fig. 2, we give the 
expected number of events per cm^ from a single GRB source, as a function 
of the cosine of zenith angle 6. 

The dependence on 6 of the muon flux is a consequence of neutrino absorption 
in Earth, which cannot be neglected at high energy. Moreover, the suppres- 
sion at cos^ < —0.8 reflects the hard discontinuity of Earth density profile 
in correspondence of the massive nucleus. In particular, in Fig. 6 the calcu- 
lated fluxes are compared with upper limits by experiments (MACRO [45], 
1MB [46], and AMANDA [29]). (Note that the 1MB limit has been lowered 
for convenience of presentation. Moreover, the AMANDA limit, shown as a 
dashed line, has been estimated starting from the declared upper limit on neu- 
trino flux). Finally, for comparison purposes, in Fig. 6 we show the number of 
events calculated by assuming as theoretical input the well-known Waxman 
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Fig. 6. The total number of events per cm^ expected from a single GRB as a function 
of the cosine of zenith angle 9. The plot is shown for the set of model parameters 
given in Fig. 2. The thin dashed line is calculated starting from the model by 
Waxman and Bahcall [25]. Upper limits by experiments MACRO [45], 1MB [46], 
and AMANDA [29] are also shown (see text for comments). 

and Bahcall ([8]) model of neutrino production from GRBs. The sensitivity 
to gamma ray bursts as point-like sources is improved by the search for space 
and time correlation with optical information. In this case the background is 
strongly suppressed. 

The same analytical procedure has been used in order to calculate the expected 
number of events induced by our predicted neutrino flux from unresolved 
GRBs labelled in Fig. 3 as GRB pN (1) and (2). The overall rate of events 
expected for E^ > 10 GeV is 190 and 46 events per km^ from the entire 
lower Earth hemisphere for the flux (1) and (2) respectively. The event rate 
for Efj^ > 10 TeV, in a range where the signal is expected to start arising 
over the atmospheric neutrino background, becomes 26 and 6 events per km^ 
respectively. In any case a good knowledge of the background is required in 
order to have a reliable measure. 

As shown in Fig. 3 and Fig. 6, current experiments have not the required 
sensitivity to exclude or confirm our predictions. Neutrinos telescopes of the 
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next generation like ANTARES [47], NEMO [48], NESTOR [49], with effective 
area ranging from 0.1 to 1 km^, will be likely capable to success in this task. 



5 Conclusions 



We have presented a model for the production of high energy neutrinos from 
GRBs in an energy range from few GeV up to few PeV, through the interac- 
tion of the shock accelerated protons with the matter surrounding the source. 
The results have been compared with other theoretical predictions and with 
the available experimental upper limits from underground and underwater/ice 
neutrino detectors. Current experiments have not the required sensitivity to 
exclude or confirm our predictions but the next generation neutrinos telescopes 
are expected to be capable to reach this goal in few years of data taking. 
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